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Abstract 

We present the numbers of spanning forests on the Sierpinski gasket SGd(n) at 
stage n with dimension d equal to two, three and four, and determine the asymptotic 
behaviors. The corresponding results on the generalized Sierpinski gasket SGd,b( n ) 
with d = 2 and 6 = 3,4 are obtained. We also derive the upper bounds of the 
asymptotic growth constants for both SGd and SG2,b- 
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1 Introduction 



The enumeration of the number of spanning forests Nsf(G) on a graph G is a 
problem of interest in mathematics [Tf2] and physics [3]. It is well known that 
the number of spanning forests is given by the Tutte polynomial T(G, x, y) 
evaluated at x — 2, y — 1 [3]. Alternatively, it corresponds to a special q — ► 
limit of the partition function of the g-state Potts model in statistical mechan- 
ics [5J. Some recent studies on the enumeration of spanning forests and the 
calculation of their asymptotic growth constants on regular lattices were car- 
ried out in Refs. [oTITf8f9][T()Ul|ll2| . It is of interest to consider spanning forests 
on self-similar fractal lattices which have scaling invariance rather than trans- 
lational invariance. Fractals are geometric structures of noninteger Hausdorff 
dimension realized by repeated construction of an elementary shape on pro- 
gressively smaller length scales [T^fTl] . A well-known example of fractal is the 
Sierpinski gasket. We shall derive the recursion relations for the numbers of 
spanning forests on the Sierpinski gasket with dimension equal to two, three 
and four, and determine the asymptotic growth constants. We shall also con- 
sider the number of spanning forests on the generalized Sierpinski gasket with 
dimension equal to two. 

2 Preliminaries 

We first recall some relevant definitions for spanning forests and the Sierpinski 
gasket in this section. A connected graph (without loops) G = (V, E) is defined 
by its vertex (site) and edge (bond) sets V and E [ToTTE] . Let v(G) = \V\ be 
the number of vertices and e(G) — \E\ the number of edges in G. A spanning 
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subgraph G' is a subgraph of G with the same vertex set V and an edge set 
E' C E. While a tree is a connected graph with no circuits, a spanning forest 
on G is a spanning subgraph of G that is a disjoint union of trees. That is, 
a subgraph of G without any cycles, or an acyclic graph. Here an isolated 
vertex is considered as a tree. The degree or coordination number ki of a 
vertex v $ G V is the number of edges attached to it. A A;- regular graph is a 
graph with the property that each of its vertices has the same degree k. In 
general, one can associate an edge weight to each edge connecting adjacent 
vertices Vi and Vj (see, for example [17]). For simplicity, all edge weights are 
set to one throughout this paper. 

When the number of spanning forests Nsf{G) grows exponentially with v(G) 
as v(G) — > oo, there exists a constant zq describing this exponential growth 



where G, when used as a subscript in this manner, implicitly refers to the 
thermodynamic limit. 

The construction of the two-dimensional Sierpinski gasket SG2{n) at stage n 
is shown in Fig.[TJ At stage n = 0, it is an equilateral triangle; while stage n+1 
is obtained by the juxtaposition of three n-stage structures. In general, the 
Sierpinski gaskets SGd can be built in any Euclidean dimension d with fractal 
dimensionality D = ln(d + 1)/ In 2 [20 J. For the Sierpinski gasket SGd(n), the 
numbers of edges and vertices are given by 
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Except the (d + 1) outmost vertices which have degree d, all other vertices of 
SGd(n) have degree 2d. In the large n limit, SGd is 2<i-regular. 



Fig. 1. The first four stages n = 0, 1,2,3 of the two-dimensional Sierpinski gasket 



The Sierpinski gasket can be generalized, denoted as SGd,b(n), by introducing 
the side length b which is an integer larger or equal to two [21J. The gener- 
alized Sierpinski gasket at stage n + 1 is constructed with b layers of stage n 
hypertetrahedrons. The two-dimensional SG2,b(n) with b = 3 at stage n = 1, 2 
and b = 4 at stage n — 1 are illustrated in Fig. [21 The ordinary Sierpinski 
gasket SGdin) corresponds to the b = 2 case, where the index b is neglected for 
simplicity. The Hausdorff dimension for SGd.b is given by D = In {^ +d d -1 J /In b 
[2T] . Notice that SGd,b is not /c-regular even in the thermodynamic limit. 

3 The number of spanning forests on SG2{n) 

In this section we derive the asymptotic growth constant for the number of 
spanning forests on the two-dimensional Sierpinski gasket SG2{n) in detail. 
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SG 2)3 {1) 



SG 2 , 3 (2) 



SG 2A (1) 



Fig. 2. The generalized two-dimensional Sierpinski gasket SG2 ) b( n ) with b = 3 at 
stage n = 1, 2 and 6 = 4 at stage n = 1. 

Let us start with the definitions of the quantities to be used. 

Definition 3.1 Consider the generalized two-dimensional Sierpinski gasket 
SG 2 ,b(n) at stage n. (a) Define f2,b( n ) = N S F(SG2,b{n)) as the number of 
spanning forests, (b) Define t2^{ n ) as the number of spanning forests such 
that the three outmost vertices belong to one tree, (c) Define ga,2,b{n), gb2,b{n), 
9C2,b{n) as the number of spanning forests such that one of the outmost vertices 
belongs to one tree and the other two outmost vertices belong to another tree, 
(d) Define h 2 ^{n) as the number of spanning forests such that each of the 
outmost vertices belongs to a different tree. 

Since we only consider ordinary Sierpinski gasket in this section, we use the 
notations f 2 (n), t 2 (n), ga 2 (n), gb 2 (n), gc 2 (n) and h2{n) for simplicity. They are 
illustrated in Fig. [31 where only the outmost vertices are shown. It is clear that 
the values ga,2{n), gb2(n), gc2{n) are the same because of rotation symmetry, 
and we define g2 {n) = ga2{n) = gb2{n) = gc2{n). Similarly for the generalized 
case, we define g 2 ,b(n) = ga2,b{n) = gb 2 ,b{p) = gc2,b(n). It follows that 



/ 2 (n) = t 2 {n) + 3g 2 (n) + h 2 {n) . 



(4) 



5 



The initial values at stage zero are £2(0) = 3, #2(0) = 1, h 2 (0) = 1 and 
^2(0) = 7. The purpose of this section is to obtain the asymptotic behavior of 
f 2 {n) as follows. The three quantities t 2 (n), g 2 (n) and h2{n) satisfy recursion 
relations. 




/ 2 (n) t 2 (n) ga 2 (n) gh(n) gc 2 (n) h 2 (n) 

Fig. 3. Illustration for the spanning subgraphs f 2 (n), t 2 (n), ga 2 (n), gb 2 (n), gc 2 (n) 
and h 2 {n). The two outmost vertices at the ends of a solid line belong to one tree, 
while the two outmost vertices at the ends of a dot line belong to separated trees. 

Lemma 3.1 For any non- negative integer n, 

f 2 (n + l) = f!(n)-lt 2 (n)+g 2 (n)f, (5) 
t 2 ( n + 1) = Gtj{n)g 2 {n) + 3t 2 {n)g 2 2 {n) , (6) 



g 2 (n + l)=t 2 2 {n)h 2 (n) + 2t 2 {n)g 2 (n)h 2 (n) + lt 2 {n)gl{n) + 4g 3 2 (n) 

+g 2 2 (n)h 2 (n) , (7) 

h 2 {n + 1) = 12t 2 {n)g 2 (n)h 2 {n) + Ug 3 2 {n) + 24g 2 2 (n)h 2 {n) + 9g 2 (n)h 2 2 {n) 

+3t 2 (n)h 2 2 (n) + h 3 2 (n) . (8) 

Proof The Sierpinski gaskets SG 2 (n + 1) is composed of three SG 2 (n) with 
three pairs of vertices identified. For the number f 2 (n + 1), the unallowable 
configurations are those with a circuit, i.e., the two identified vertices of each 
SG 2 (n) belong to the same tree as illustrated in Fig. HI Therefore, we have 

f 2 (n + 1) = /|(n) - [t 2 (n) + ga 2 {n)][t 2 {n) + gb 2 {n)][t 2 {n) + gc 2 {n)] . (9) 
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With the identity ga2(n) = gb2(n) = gc2(n) = (^(ft), Eq. (|3j) is verified. 




"ZV;,-,:XA",:XZ;, 

Fig. 4. Illustration for the expression of f2(n + 1)- 

The number ^(ft + 1) consists of six configurations where two of the SG2(n) 
are in the t2(n) status and the other one is in the (n) status, and three 
configurations where one of the SG2{n) is in the £2 (ft) status and the other 
two are in the 52 (ft) status as illustrated in Fig. Therefore, we have 

t 2 (n + l) = 2tl(n)[ga 2 {n) + gb 2 (n) + gc 2 {n)} + t 2 {n)ga 2 {n)gb 2 {n) 

+t 2 {n)ga 2 (n)gc 2 {n) +t 2 {n)gb 2 {n)gc 2 {n) . (10) 

With the identity ga 2 {n) = gb 2 {n) = gc2(n) = #2 (ft), Eq. (Q is verified. 




Fig. 5. Illustration for the expression of t 2 (n + 1). The multiplication of three on 
the right-hand-side corresponds to the three possible orientations of SG 2 (n + 1). 

Similarly, ga 2 (n + 1) for SG 2 (n + 1) can be obtained with appropriated con- 
figurations of its three constituting SG 2 (n) as illustrated in Fig. El Thus, 

ga 2 (n + 1) =tl(n)h 2 (n) + t2(n)ga 2 (n)[ga 2 (n) + gc 2 {n) + h 2 (n)} 

+t 2 {n)ga2{n)[ga 2 {n) + gb 2 (n) + h 2 (n)} + f2{n)ga 2 2 (n) 
+t 2 {n)ga 2 {n)gc 2 {n) + t 2 (n)ga 2 (n)gb 2 (n) 



7 



+ga 2 (n)gb 2 (n)gc 2 (n) . (11) 

With the identity ga2{n) = gbzin) = gc2(n) = g2(n) and Eq. (J4j), Eq. (JTj) is 
verified. 





+ 



+ 




Fig. 6. Illustration for the expression of ga 2 (n + 1). 



Finally, h 2 (n+l) is the summation of appropriated configurations as illustrated 
in Fig. [3, so that 

h 2 {n + 1) 

= 4:t 2 (n)h 2 (n)[ga 2 (n) + gb 2 (n) + gc 2 {n)} 

+2gc 2 (n)ga 2 (n)[gc 2 (n) + ga 2 (n)] + 2ga 2 (n)gb 2 (n)[ga 2 {n) + gb 2 (n)} 
+2gb 2 (n)gc 2 (n)[gb 2 (n) + gc 2 (n)} + 2ga 2 (n)gb 2 (n)gc 2 (n) 
+3t 2 (n)hl(n) + ft|(n) + 3[ga 2 (n) + ^(n) + gc 2 (n)]h\(n) 
+{3[ga 2 (n) + gb 2 (n) + gc 2 (n)] 2 - ga 2 2 {n) - gb 2 2 {n) - gc 2 2 {n)}h 2 {n) .(12) 

With the identity ga 2 (n) = gb 2 {n) = gc 2 {n) = g 2 (n), Eq. (jSj) is verified. 

Eq. (ED can also be obtained by substituting Eqs. (jSJ), (CJ) and (jSJ) into Eq. (@J. 

□ 

The values of f 2 (n), t 2 (n), g 2 (n), h 2 {n) for small n can be evaluated recursively 
by Eqs. (j5J), (EI), (J7J), (jHJ) as listed in Tabled! These numbers grow exponentially, 
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Fig. 7. Illustration for the expression of h2(n + 1). The multiplication of three on 
the right-hand-side corresponds to the three possible orientations SG2(n + 1). 
and do not have simple integer factorizations, in contrast to the corresponding 

results for the number of spanning trees [22J. To estimate the value of the 
asymptotic growth constant defined in Eq. ((T]), we need the following lemma. 

Table 1 

The hrst few values of fain), t2(n), g2(n), h^in). 



n 





1 


2 


3 


h(n) 


7 


279 


20,592,775 


8,696,126,758,781,951,722,199 


t 2 (n) 


3 


63 


1,294,083 


36,212,372,367,917,382,063 




1 


41 


2,022,893 


215,741,040,104,979,715,185 


h 2 (n) 


1 


93 


13,230,013 


8,012,691,266,099,095,194,581 



Lemma 3.2 The asymptotic growth constant for the number of spanning 
forests on 5 , G 2 (n) is bounded: 



2 2 
— ln/i 2 (m) < z SG2 < _lB/ 2 (m) 



(13) 



where m is a positive integer. 
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Proof We first show that the ratio t2(n)/g2(n) is a strictly decreasing se- 
quence. By Eqs. (jHj) and (El), we have 



t 2 {n + l) 
92(n + 1) 

6t 2 2 {n)g 2 {n) + 3t 2 (n)g 2 (n) 



t|(n)/i 2 (n) + 2t 2 (n)^ 2 (n)/i 2 (n) + 7t 2 (n)g 2 (n) + 4(/f(n) + ^ 2 (n)/i 2 (n) 
t 2 (n)[Gt 2 (n)g 2 (n)+3g 2 (n)] Gt 2 (n) 

g2(n)[7t2(n)g 2 (n)+4gl(n)\ 7g 2 (n) ' { ' 



From the values in Table (TJ t 2 (n) / g 2 (n) is less than one for n > 1. It is clear 
that this ratio approaches to zero as n increases. Similarly, g 2 (n)/h 2 (n) is also 
a strictly decreasing sequence by Eqs. ([7j) and (jSJ). 



fl 2 (rc + l) 
/i 2 (n + 1) 



< 



< 



3t 2 (w)g 2 (n)fr 2 (n) + 7t 2 (n)^ 2 2 (n) + 4gl(n) + gj(n)h 2 (n) 

12t 2 (r% 2 (n)/i 2 (n) + Ugi(n) + 2Ag 2 {n)h 2 {n) + 9g 2 (n)h 2 2 (n) + 4t 2 (n)h 2 2 (n) 

92(n)[St 2 (n)h 2 (n) + 7t 2 (n)g 2 (n) + Ag 2 (n) + g 2 (n)h 2 (n)] 
h 2 (n)[it 2 (n)h 2 (n) + 12t 2 (n)g 2 (n) + 24g 2 (n) + 9g 2 (n)h 2 (n)} 

4/i 2 (n) 

where we have used the fact that t 2 (n) < g 2 (n) < /i 2 (n) for n > 1. Again, 
g 2 (n)/h 2 (n) approaches to zero as n increases. The relation t 2 (n) <C g 2 {n) <C 
/t 2 (n) for large n is expectable since it is rare to keep the three outmost 
vertices of SG 2 (n) in the same tree for t 2 (n) and h 2 (n) should dominate when 
n becomes large. In fact, both f 2 (n) and g 2 {n) are negligible compared with 
h 2 (n) such that f 2 (n) ~ h 2 (n) for large n. By Eqs. (jSJ) and (jSJ), we have the 
upper and lower bounds for f 2 (n): 



h\{n - 1) < fr 2 (n) < / 2 (n) < /f(n - 1) , (16) 



10 



such that 

h 2 {mf n - m < f 2 (n) < / 2 (m) 3 "" m , (17) 

where m is a fixed integer. With the definition for z$g 2 given in Eq. ([I]) and the 
number of vertices of SG2{n) is 3(3™ + l)/2 by Eq. ([3]), the proof is completed. 
□ 

Proposition 3.3 The asymptotic growth constant for the number of span- 
ning forests on the two-dimensional Sierpinski gasket SG2{n) in the large n 
limit is z SG2 = 1.24733719931.... 

Proof Define ratios a{n) = t 2 (n) / / 2 (n) and (3{n) = g 2 (7i)// 2 (n). By Eq. (jlj), 



it is clear that < a(n) +(3(n) < 1. According to Lemma |3~2| a(n) + (3(n) is a 
strictly decreasing sequence. By Eq. ([5]), let us define r{n) = f 2 {n) / f^in — l) = 
1 — [a(n — 1) + f3(n — l)] 3 for positive integer n . It follows that 

ln/2(n) =31n/ 2 (7i — 1) + lnr(n) = ... 

n 

= 3 n - m ln/ 2 (m) + Y, 3 n ~ j In r(j) 

j=m+l 

(on—m 1 \ 
jlnr(m + l). (18) 

Divide this equation by 3(3™ + l)/2 and take the limit n — > 00, the difference 
between the upper bound in Eq. (JIB"]) and the asymptotic growth constant is 
bounded: 

^ln/ 2 (m)-^ G2 <^ T ln(l-[«(m)+/3(m)] 3 ) . (19) 

When m is as small as three, the right-hand-side of Eq. (jT9l) is about 3 x 10~ 7 
by the values given in Table [TJ Similarly, it can be shown that the difference 
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between zsg 2 an d the lower bound (left-hand-side of Eq. ([13]) ) quickly con- 
verges to zero as m increases. In another word, the numerical values of In /2(w) 
and ln/i 2 (^) are almost the same except the first few m, and the upper and 
lower bounds in Eq. (Jl3"j) converge to the quoted value of zsg 2 - m fact, one ob- 
tains the numerical value of z$g 2 with more than a hundred significant figures 
accurate when m is equal to eight. □ 

4 The number of spanning forests on SG2,b(n) with b = 3,4 

The method given in the previous section can be applied to the number of 
spanning forests on SGd,b{ n ) with larger values of d and b. The number of 
configurations to be considered increases as d and b increase, and the recursion 
relations must be derived individually for each d and b. In this section, we 
consider the generalized two-dimensional Sierpinski gasket SG2,b{ n ) with the 
number of layers b equal to three and four. For SG2^(n), the numbers of edges 
and vertices are given by 



where the three outmost vertices have degree two. There are (6 n — 1)/5 vertices 
of SG2,3{n) with degree six and 6(6™ — l)/5 vertices with degree four. By 
Definition l3.lt the number of spanning forests is f2,z{ n ) — ^2,3 i n ) + 3<?2,3 i n ) + 
/i 2)3 (n). The initial values are the same as for SG2: £2,3(0) = 3, #2,3(0) = 1, 
^2,3(0) = 1 and ^2,3(0) = 7. By the method illustrated in the previous section, 
we obtain following recursion relations for any non-negative integer n. 



e(SG 2>3 (n)) = 3x6 



(20) 



v(SG 2 , 3 (n)) 



7 x 6 n + 8 
5 



(21) 
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/2, 3 (n + l) 

= /!» - ZfUn)[t^{n) + g 2 , 3 (n)] 3 - 3f 2 , 3 (n)[t 2 , 3 (n) + g 2 , 3 (n)] 5 

-[hA n ) + 02,3 Wf + 6t 2 ,3M/2 1 3M[*2 1 3(rO + 02,3 W] 4 

+6*l )3 (n)[* 2)3 (n) + g 2 , 3 (n)] 4 - 6t| 3 (n)[t 2 , 3 (n) + ^, 3 (n)] 3 , (22) 
*2 )3 (n + l) 

= 142^ )3 (n) ff | >3 (n) + 18t 4 3 (n)(? 2i3 (n)/i 2i3 (n) + 45^ 3 (n)4 3 (n)/i 2i3 (n) 
+153i| 3 (n)4 3 (n) + 36^ )3 (n)4 3 (n)/i 2 , 3 (n) + A5t 2 , 3 (n)gl 3 (n) 
+^{n)g{ 3 {n)h 2t3 {n) + 2^ 2 6 3 (n) , (23) 

52,3(^ + 1) 

= 77tl 3 (n)gl 3 (n)h 2 , 3 (n) + 171t| 3 (n)(7 2 4 3 (n) + 2t 4 3 (n)/4 3 (n) 
+22^ i3 (n)^ 2 , 3 (n)/i| 3 (n) + 200t| 3 (n)(? 2 3 3 (n)/i 2i3 (n) + 195t 2i3 (n)^ )3 (n) 
+4,3( n ) h h( n ) + 50^ i3 (n)4 3 (n)/i| 3 (n) + 169t 2i3 (n)^n)/i 2 , 3 (n) 
+56^ >3 (n) + 3tl i3 (n)^ 2)3 (n)^ ) 3 (n) + 42* 2)3 (n)^ i3 (n)^ i3 (n) 
+46<? 2 5 3 (n)/i 2 , 3 + 3t 2 , 3 (n)4 3 (n)/i 3 3 (n) + 12^ 4 3 (n)^ 3 (n) 
+02>)^2>) , (24) 

/i 2 ,3(n + l) 

= 60^ i3 (n)^ 2)3 (n)/i| 3 (n) + 564t| 3 (n)4 3 (n)/i 2;3 (n) + 468t 2i3 (n) ff | )3 (n) 
+14t 3 3 (n)/i 3 3 (n) + 552^ 3 (n)4 3 (n)/i| 3 (n) + 1608t 2i3 (n)4 3 (n)/i 2)3 (n) 
+468«? 2 6 3 (n) + 162t| 3 (n)( ?2 , 3 (n)/i 3 3 (n) + 1404t 2 , 3 (n)4 3 (n)^ 3 (n) 
+12364 3 (n)/i 2)3 (n) + 15^ ;3 (n)/i 4 3 (ra) + 522t 2 , 3 (n)3 2 2 3 (n)/4 3 (n) 
+11524 3 (n)^ 3 (n) + m 2t3 {n)g 2 , 3 {n)h\^n) + 534<? 2 3 3 (n)^ >3 (n) 
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+6t 2>3 (n)h 5 2)3 (n) + 135^ i3 (n)/4 )3 (n) + lSg 2 , 3 (n)h 5 2)3 (n) + h 6 2>3 (n) . (25) 

The figures for these configurations are too many to be shown here. Some 
values of f2,3(n), £2,3 (^), #2,3 (n), ^2,3 (n) are listed in Table [2j These numbers 
grow exponentially, and do not have simple integer factorizations. 
Table 2 

The first few values of f2,3{n), t 2i3 (n), ^2,3(^)1 h 2 ^{n). 



n 





1 


2 


/2,3(X) 


7 


61,905 


53,145,523,900,850,102,434,114,604,001 


*2,3(n) 


3 


8,372 


218,891,276,004,139,532,538,695,680 


52,3 (n) 


1 


8,020 


1,242,664,072,161,818,527,545,741,824 


7*2,3 (ji) 


1 


29,473 


49,198,640,408,360,507,318,938,682,849 



By a similar argument as Lemma [3T2l the asymptotic growth constant for the 
number of spanning forests on SG 2< 3(n) is bounded: 

5 5 
77"^ ln M™) < z sg*,s < In / 2 ,3(m) , (26) 

with m a positive integer. We have the following proposition. 

Proposition 4.1 The asymptotic growth constant for the number of span- 
ning forests on the two-dimensional Sierpinski gasket SG2,s{n) in the large n 
limit is z S g 2<3 = 1.31235755933.... 

The convergence of the upper and lower bounds remains quick. By the same 
method as given in the proof of Proposition 13.31 the difference between the 
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upper bound in Eq. (|26|) and the asymptotic growth constant is bounded: 



7 x 6 m 



5 



ln/ 2i3 (m) - z sch < 



7x6 



-1 



rn 



In 1-7 




(27) 



More than a hundred significant figures for zsg 2 , 3 can be obtained when m is 



equal to five. 



For S , G 2i4 (n), the numbers of edges and vertices are given by 



e(SG 2A {n)) = 3 x 10 



(28) 



v(SG 2A (n)) 



4 x 10 n + 5 
3 



(29) 



where again the three outmost vertices have degree two. There are (10 n — l)/3 
vertices of SG 2t 4:(n) with degree six, and (10 n — 1) vertices with degree four. By 
Definition 13. 1[ the number of spanning forests is / 2)4 (n) = t 2 , 4 (^) + 3g 2i4 (n) + 
/i 2i4 (n). The initial values are the same as for SG2: ^2,4(0) = 3, #2,4(0) = 1, 
^2,4(0) = 1 and ^2,4(0) = 7. We write a computer program to obtain the 
recursion relations. Using the shorthand notation tg2^{n) = t2^{n) + (72,4(^)5 
we have 



/ 2 ,4(n + l) 

= fl%n) - QflMtgUn) - 9fl,(n)tgl 4 (n) + l%f 2A {n)tg 4 2A {n)t 2A {n 
+yt A {n)tgl A {n) + 18/ 2 4 4 (n)^ 2 4 4 (n)t| 4 (n) - 18/ 2 4 4 (n)^ 2 3 4 (n)4 4 (n; 
-Qf! A {n)tg 7 2A {n) + 30/ 2 3 4 (n)^ 2 6 4 (n)t 2j4 (n) - 30/ 2 3 4 (n)^ 2 5 4 (n)t 2 A ( n ] 
+3/ 2 2 4 (n)^ 2 8 4 (n) + 2Afl A {n)tgl A {n)t2 , 4 (n) - 36/ 2 2 4 (n)^ 2 6 4 (n)t| 4 (n; 
-hAfl A {n)tgl A {n)tl 4 (n) + 60/ 2 2 4 (n)^ 2 4 4 (n)t 4 4 (n) - 5/ 2 , 4 (n)^ 2 9 4 (n; 
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+A2f 2A (n)tg 8 2A (n)t 2A (n) - A2f 2A (n)tg 7 2A (n)t 2 2A (n) 
-imf 2A {n)tg 6 2A {n)t 3 2A {n) + 330 f 2A (n)tg 5 2A (n)4 A (n) 
-162 f 2A {n)tgl A {n)t 5 2A {n) + 8f 2>4 (n)tg 3 2A (n)t 6 2A (n) + A2tg 8 2A (n)t 2 2A (n) 
-162tg 7 2A {n)t 3 2A {n) + 102tg 6 2A {n)4 A {n) + 288tgl A (n)t\ A (n) 
-A32tgt A {n)tlM) + lQ2tg 3 A (n)t 7 A (n) . (30) 

The other recursion relations for SG 2A (n) are too lengthy to be included here. 
They are available from the authors on request. Some values of f 2 ^{n), t 2 ^{n) : 
g 2 ^(n), h 2 ^{n) are listed in Table [31 These numbers grow exponentially, and 
do not have simple integer factorizations. 



Table 3 

The first few values of f 2 A(n), t 2A (n), #2,4(^)1 h 2 ^(n). 



n 


1 


2 


f 2,4(n) 


75,908,209 


6,053,025,303,996,636,848,970,430,785,675,468,144,409,657,412,247,800,423,390,303,465,602,821,564,523,873 


t 2A (n) 


6,665,475 


772,069,425,849,585,011,183,346,692,712,538,703,294,972,628,973,372,161,275,424,155,207,555,217,357 




8,406,453 


17,447,838,129,920,655,302,865,270,986,884,479,355,572,603,291,172,150,410,900,983,156,421,717,259,395 


h 2A (n) 


44,023,375 


5,999,909,720,181,025,298,050,651,626,022,102,167,639,644,629,745,310,599,996,325,091,978,348,857,528,331 



By a similar argument as Lemma [3T2l the asymptotic growth constant for the 
number of spanning forests on SG 2A (n) is bounded: 

3 3 

\nh 2A (m) < z S G 2i < -. ——\nf 2A (m) , (31) 



4 x 10 m ' v ' 4 x 10 r 



with m a positive integer. We have the following proposition. 
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Proposition 4.2 The asymptotic growth constant for the number of span- 
ning forests on the two-dimensional Sierpinski gasket SG2,4(n) in the large n 
limit is z SG2A = 1.36051646575.... 



The convergence of the upper and lower bounds is again quick. By the same 
method as given in the proof of Proposition 13. 3[ the difference between the 
upper bound in Eq. (13~TT) and the asymptotic growth constant is bounded: 



4 x 10 r 



ln/ 2l4 (m) - z SG2A < 



-1 



12 x 10 r 



In 1 - 15 



tg2A m ) 



(32) 



More than a hundred significant figures for zsg 2 a can be obtained when m is 
equal to four. 



5 The number of spanning forests on SGd{n) with d — 3, 4 

In this section, we derive the asymptotic growth constant of spanning forests 
on SGd{n) with d = 3,4. For the three-dimensional Sierpinski gasket SG^n), 
we use the following definitions. 

Definition 5.1 Consider the three-dimensional Sierpinski gasket SG^(n) at 
stage n. (a) Define fz{n) = NsFiSG^^n)) as the number of spanning forests, 
(b) Define £3(71.) as the number of spanning forests such that the four outmost 
vertices belong to one tree, (c) Define (n) as the number of spanning forests 
such that one of the outmost vertices belongs to one tree and the other three 
outmost vertices belong to another tree, (d) Define h^n) as the number of 
spanning forests such that two of the outmost vertices belong to one tree and 
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the other two outmost vertices belong to another tree, (e) Define p 3 (n) as the 
number of spanning forests such that two of the outmost vertices belong to one 
tree and the other two outmost vertices separately belong to other trees, (f) 
Define q 3 (n) as the number of spanning forests such that each of the outmost 
vertices belongs to a different tree. 

The quantities f 3 (n), t 3 (n), g 3 (n), h 3 (n), p 3 (n) and (73 (n) are illustrated in 
Fig. [H], where only the outmost vertices are shown. There are four equivalent 
#3 (n), three equivalent h 3 (n), and six equivalent p 3 {n). By definition, 

f 3 {n) = t 3 {n) + 4g 3 (n) + 3h 3 {n) + 6p 3 (n) + q 3 {n) . (33) 

The initial values at stage zero are £3(0) = 16, #3(0) = 3, h 3 (0) = 1, ^3(0) = 1, 
g 3 (0) = 1 and / 3 (0) = 38. 




/3M t 3 (n) g 3 (n) h 3 (n) p 3 (n) q 3 (n) 

Fig. 8. Illustration for the spanning subgraphs f 3 (n), t 3 (n), g 3 (n), h 3 (n), p 3 (n) and 
q 3 (n). The two outmost vertices at the ends of a solid line belong to one tree, while 
the two outmost vertices at the ends of a dot line belong to separated trees. 

The recursion relations are lengthy and given in the appendix. Some values 
of f 3 {n), t 3 (n), g 3 (n), h 3 (n), p 3 (n), q 3 {n) are listed in Table HI These numbers 
grow exponentially, and do not have simple integer factorizations. 
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Table 4 

The first few values of fo(n), t 3 (ri), g 3 {n), h 3 (n), pz(n), qz{n). 



n 





1 


2 


h(n) 


38 


701,866 


150,308,440,552,729,541,599,408 


Un) 


16 


173,880 


14,568,001,216,879,127,537,520 


93(11) 


3 


63,354 


10,109,099,387,983,187,560,398 


h 3 (n) 


1 


9,059 


1,150,970,295,799,746,536,513 


ps(n) 


1 


31,357 


9,282,357,698,529,097,198,747 


93(«) 


1 


59,251 


36,156,984,705,343,841,018,275 



By a similar argument as Lemma I3~2l the asymptotic growth constant for the 
number of spanning forests on SGs(n) is bounded: 



In q 3 (m) ln/ 3 (m) 
2 x 4m ZsG * 2x4™ 



(34) 



with m a positive integer. We have the following proposition. 



Proposition 5.1 The asymptotic growth constant for the number of span- 
ning forests on the three-dimensional Sierpinski gasket SG 3 (n) in the large n 
limit is z SG3 = 1.66680628117.... 



The convergence of the upper and lower bounds is not as quick as for the two 
dimensional cases. By the same method as given in the proof of Proposition 
13. 3[ the difference between the upper bound in Eq. (1341) and the asymptotic 
growth constant is bounded: 
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1 



2x4 m 



ln/ 3 (m) - z SG3 



< 



1 



6 x 4 m 



In 1-7 



i 3 (m) 2g 3 (m) fe 3 (m) p 3 (m) 



1 3^ 



/3H /s("l) /3("l) M™) 



(35) 



More than a hundred significant figures for z S g 3 can be obtained when m is 
equal to nine. 



For the four- dimensional Sierpinski gasket SG±{n), we use the following defi- 
nitions. 



Definition 5.2 Consider the four-dimensional Sierpinski gasket SG^(n) at 
stage n. (a) Define / 4 (n) = NsFiSG^n)) as the number of spanning forests, 
(b) Define £ 4 (n) as the number of spanning forests such that the five outmost 
vertices belong to one tree, (c) Define g 4 (n) as the number of spanning forests 
such that two of the outmost vertices belong to one tree and the other three 
outmost vertices belong to another tree, (d) Define /i 4 (n) as the number of 
spanning forests such that one of the outmost vertices belong to one tree and 
the other four outmost vertices belong to another tree, (e) Define p 4 (n) as the 
number of spanning forests such that one of the outmost vertices belong to 
one tree, two of the other outmost vertices belong to another tree and the rest 
two outmost vertices belong to a third tree, (f) Define g 4 (n) as the number 
of spanning forests such that three of the outmost vertices belong to one tree 
and the other two outmost vertices separately belong to other trees, (g) Define 
r 4 (n) as the number of spanning forests such that two of the outmost vertices 
belong to one tree and the other three outmost vertices separately belong to 
other trees, (h) Define s 4 (n) as the number of spanning forests such that each 
of the outmost vertices belongs to a different tree. 
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The quantities f 4 (n), t 4 (n), ^(n), h 4 (n), p 4 (n), q±(n), r 4 (n) and s 4 (n) are 
illustrated in Fig. M, where only the outmost vertices are shown. There are ten 
equivalent g 4 (n), five equivalent h 4 (n), fifteen equivalent p 4 {n), ten equivalent 
q 4 (n) and ten equivalent r 4 (n). By definition, 

f A (n) = t 4 (n) + 10g 4 (n) + 5/i 4 (n) + 15p 4 (n) + 10g 4 (n) + 10r 4 (n) + s 4 (n) . 

(36) 

The initial values at stage zero are t 4 (0) = 125, g 4 (0) = 3, h 4 (0) = 16, p 4 (0) = 
1, g 4 (0) = 3, r 4 (0) = 1, s 4 (0) = 1 and / 4 (0) = 291. 




h{n) t 4 (n) g 4 {n) h 4 (n) p 4 {n) 




q 4 (n) r 4 (n) s 4 (n) 

Fig. 9. Illustration for the spanning subgraphs / 4 (n), t 4 (n), g 4 (n), h 4 (n), p 4 (n), 
q 4 (n), r 4 {n) and s 4 (n). The two outmost vertices at the ends of a solid line belong 
to one tree, while the two outmost vertices at the ends of a dot line belong to 
separated trees. 

We write a computer program to obtain the recursion relations. Using the 
shorthand notations tr 4 {n) = t 4 {n) + 4g 4 (n) + 3h 4 (n) + 3p 4 (n) +3g 4 (n) +r 4 (n), 
tq 4 (n) = t 4 (n)+g 4 (n)+2h 4 (n) + q 4 (ri), tp 4 (n) = t 4 (n) +2g 4 (n) + h 4 (n) +p 4 (n) 
and th 4 {n) = t 4 {n) + h 4 (n), we have 

f 4 (n + l) 

= f!(n) - 10 /f(n)tr|(n) - 15f 4 {n)tr 4 4 {n) - 30/ 4 (n)tg 4 4 (n) - 12tr 4 5 (n) 
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+mf 4 (n)trl(n)tql(n) - 15tri(n)tp 4 (n) + 30trj(n)th 4 (n) 
-30tr'i(n)tpi(n) + 120tr'i(n)tp 4 (n)th 4 (n) + UOtrl(n)tqj(n) 
-120trl(n)thl(n) + 2A0trl(n)tql(n)tp 4 (n) - 480trl(n)tqj(n)th 4 (n) 
+300tr 4 (n)tql(n)tp 2 4 (n) - 1200tr 4 (n)tql(n)tp 4 (n)th 4 (n) 
-180tn(n)tqt(n) + I2mtr 4 (n)tql(n)th\(n) - 51tpf(n) 
+510tpi(n)th 4 (n) + 2mtql(n)tp\(n) - 2040t#[(n)t/i 4 (n) 
-1560tql{n)tpl(n)th 4 (n) + 4080tpl{n)thl{n) - 210tg|(n)tp 4 (n) 
+3120tql(n)tp 4 (n)thl(n) - A080tp 4 {n)th\{n) + A20tql(n)th 4 {n) 
-2080tql(n)thl(n) + 1632th 5 4 {n) . (37) 

The other recursion relations for SG 4 (n) are too lengthy to be included here. 
They are available from the authors on request. Some values of f 4 (n), t 4 (n), 
g 4 {n), h 4 (n), p 4 (n), q 4 (n), r 4 (n), s 4 (n) are listed in Table [5j These numbers 
grow exponentially, and do not have simple integer factorizations. 



By a similar argument as Lemma [3T2l the asymptotic growth constant for the 
number of spanning forests on SG 4 (n) is bounded: 

2 2 

■ In s 4 (m) < z SGa < — — T ^ f 4 (m) , (38) 



with m a positive integer. We have the following proposition. 



Proposition 5.2 The asymptotic growth constant for the number of span- 
ning forests on the two-dimensional Sierpinski gasket SG 4 (n) in the large n 
limit is z SGi = 1.98101707560.... 
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Table 5 

The first few values of f^n), tn(n), 54 (n), h^n), pi(n), q^n), r^n), 54 (n). 



n 


1 


2 


kin) 


85,824,132,029 


7,035,17,527,028,105,500,700,677,412,563,863,619,648,991,055,157,831,483 




3,412,986,435 


96,263,552,482,319,683,899,326,687,304,651,572,426,360,843,549,870,965 


94 (») 


392,122,089 


2,066,883,222,491,708,347,294,489,449,954,683,350,540,164,424,914,435 


hi{n) 


5,923,774,096 


40,841,537,587,690,687,322,887,835,686,137,425,636,710,177,922,212,520 


p 4 (n) 


224,652,411 


1,952,486,255,633,069,494,764,677,365,066,319,434,639,193,908,980,317 


g 4 (n) 


1,740,690,487 


19,621,800,909,697,266,778,177,200,667,594,598,639,201,513,851,821,683 


r4(n) 


693,438,141 


12,210,477,454,458,190,580,945,663,798,559,596,025,810,422,699,074,029 


s 4 (n) 


1,159,981,779 


34,767,376,906,364,680,701,267,847,191,441,347,363,970,403,604,091,693 



The convergence of the upper and lower bounds is even slower compared with 
that for SG3. By the same method as given in the proof of Proposition [33J the 
difference between the upper bound in Eq. ( 13"5j) and the asymptotic growth 
constant is bounded: 



ln/ 4 (m) - z SGi < 



-1 
2 x 5 m+1 



In 1 - 67 



£r 4 (m) 



(39) 



We only have fourteen significant figures for zsg^ with m calculated up to six. 
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6 Bounds of the asymptotic growth constants 



As the spanning tree is a special case of spanning forest where there is only 
one component, it is clear that the number of spanning trees Nst{G) is always 
less than Nsf{G). Define 

z G = hm , 40 

v{G)-*oo V(&) 

then z G < zq. We have obtained such asymptotic growth constants for the 
number of spanning trees on the Sierpinski gasket SGd for general d and SG2,b 
with b = 3,4 in Ref. [22]. They serve as the lower bounds for our current 
consideration for the spanning forests. 



By Eq. ([3]) and a similar argument as Lemma 13721 we have the upper bound of 
the asymptotic growth constant for the number of spanning forests on SGd(n): 

ZSG * < ( d + \)m + i lnN SF{SG d {m)) , (41) 

with m a positive integer. Although the number NsF{SGd{m)) for general m 
is difficult to obtain, it is known for m = 0. We first recall that SGd{0) at stage 
zero is a complete graph with (d + 1) vertices, each of which is adjacent to all 
of the other vertices. The number of spanning forests on the complete graph is 
given by sequence A001858 in Ref. [23]. The first few values of NsF(SGd(0)) 
are 7, 38, 291, 2932 for d from 2 to 5 1241. Define 



zsG d = j^\nN SF (SGd(0)) , (42) 



then zsc d < z sc d - We list the first few values of z SGd , zsG d , z SG d an d their 
ratios in Table El Notice that the upper bound is closer to the exact value 
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Table 6 



Numerical values of z SGd , zsG d , z SG d an d their ratios. The last digits given are 
rounded off. 



d 


D 


Z-SG d 


z SG d 


*SG d 


lSG d l z SG d 


ZSG d /zsG d 


2 


1.585 


1.048594857 


1.247337199 


1.297273433 


0.8406667076 


0.9615067787 


3 


2 


1.569396409 


1.666806281 


1.818793080 


0.9415589724 


0.9164353546 


4 


2.322 


1.914853265 


1.981017076 


2.269329307 


0.9666010902 


0.8729526691 


5 


2.585 


2.172764568 




2.661146688 







when d is small, while the lower bound is closer to the exact value when d is 
large. 



For the generalized Sierpinski gasket SG2j>{ n ) with dimension equal to two, 
the number of vertices can be calculated to be 



v(SG 2tb (n)) 



6 + 4 
6 + 2 



6(6+1) 



2(6+1) 
6 + 2 



(43) 



The upper bound of the asymptotic growth constant for the number of span 
ning forests on SG2,b(n) is given by 

'6 + 2\ ]nN SF {SG 2 , b (m)) 



ZSG 2 , b < 



b + 4j [6(6 + l)/2]' 



(44) 



with m a positive integer. Although the number NsF{SG2,b{fn)) for general m 
is difficult to obtain, it is always equal to seven for stage zero since SG 2t b(^) 
is the equilateral triangle. Define 

6 + 2 



Z.SG: 



2,6 



6 + 4 



In 7 , 



(45) 
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Table 7 



Numerical values of z_sg 2 , zsg 2 &j s&Ga 6 an d their ratios. The last digits given are 
rounded off. 



b 


D 




z SG 2fi 


z SG 2>b 


z SG 2 J z SG 2 , h 


ZSG 2 , b /zSG 2 , b 


3 


1.631 


1.133231895 


1.312357559 


1.389935821 


0.8635084908 


0.9441857241 


4 


1.661 


1.194401491 


1.360516466 


1.459432612 


0.8779030028 


0.9322228754 


oo 


2 






1.945910149 







then z S G 2 , b < ZSG 2 , b - We list the first few values of z SG2 b , z SG2 b , z SG%b and 
their ratios in Table Notice that the upper bound is closer to the exact 
value when b is small, while the lower bound is closer to the exact value when 
b is large. 



A Recursion relations for SG 3 (n) 



We give the recursion relations for the three-dimensional Sierpinski gasket 
SG^(n) here. For any non- negative integer n, we have 

h(n + l) 

= / 3 4 (n) - 4/ 3 (n)[t 3 (n) + 2g 3 (n) + h 3 (n) + p 3 (n)] 3 
-3[t 3 (n) + 2g 3 (n) + h 3 {n) + p 3 {n)f 
+12[t 3 (n) + 2g 3 {n) + h 3 {n) + p 3 {n)] 2 [t 3 {n) + g 3 (n)f 
-6[t 3 (n) + g 3 {n)}\ (A.l) 



t 3 (n+l) 
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= 72t 2 3 {n)p 3 {n)[g 3 {n) + h 3 (n)} + 56t 3 (n)[g 3 (n) + h 3 (n)} 3 + 2At 2 3 (n)p 2 (n) 
+12t 3 (n)p 3 (n)[llg 2 3 (n) + 12g 3 (n)h 3 (n) + h 2 3 {n)\ 
+I2g 2 3 {n)[3g 2 3 {n) + Sg 3 (n)h 3 (n) + 6h 2 3 (n)} 
+12t 3 (n)p 2 3 (n)[4g 3 (n) + h 3 {n)} + A8g 2 3 (n)p 3 (n)[g 3 (n) + h 3 {n)] 
+At 3 (n)p 3 3 (n) + I2g 2 3 {n)p 2 3 {n) , (A.2) 

9s(n + 1) 

= 6t 2 3 (n)q 3 (n)[g 3 (n) + h 3 {n)\ + 24t 2 3 (n)p 2 3 (n) 
+12t 3 (n)p 3 (n)[9g 2 3 (n) + 16g 3 (n)h 3 (n) + 7h 2 3 (n)} 
+4g 3 {n)[5gl{n) + 18g 2 3 (n)h 3 (n) + 24g 3 (n)h 2 3 (n) + Uh 3 3 (n)} 
+Qt 2 3 (n)p 3 (n)q 3 (n) + Qt 3 (n)p 2 3 (n)[21g 3 (n) + llh 3 (n)\ 
+3t 3 {n)q 3 {n)[5g 2 3 {n) + 6g 3 (n)h 3 (n) + h 2 3 {n)} 
+2Ag 3 (n)p 3 (n)[5gj(n) + 9g 3 {n)h 3 {n) + 4h 2 3 (n)] 
+6t 3 (n)p 3 (n)q 3 (n)[3g 3 (n) + h 3 (n)} + 2lt 3 {n)p\{n) 
+6g 3 (n)p 2 3 (n)[19g 3 (n) + 13h 3 (n)\ 

+3g 3 (n)q 3 (n)[3gl(n) + Ag 3 (n)h 3 (n) + h\{n)\ + 3t 3 (n)p 2 3 {n)q 3 (n) 

+6g 3 (n)p 3 (n)q 3 (n)[2g 3 (n) + /i 3 (n)] + 2hg 3 (n)p\{n) 

+3g 3 (n)p 2 3 (n)q 3 (n) , (A.3) 

/i 3 (n + l) 

= 2t 2 3 (n)p 2 3 (n) + 12t3(n)p 3 (^)[^(^) + 4<7 3 (n)/i 3 (n) + 3h 2 3 (n)\ 
+2[g 4 3 (n) + 8«? 3 3 (n)/i 3 (n) + 18g 2 3 (n)h 2 3 (n) + Wg 3 (n)h 3 3 (n) + 11/if (ra)] 
+8* 3 (ra)^(n)[20 3 (n) + 3/i 3 (ra)] 

+8p 3 (^)[2(?3 3 (n) + 9^(n)/i 3 M + 9^(n)/i 2 (n) + 2/i 3 (n)] 
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+U 3 {n)p\{n) + 2pl(n)[10gl(n) + 24g 3 (n)h 3 (n) + 9h 2 3 {n)\ 

+8p 3 3 (n) [g 3 (n) + h 3 (n)\ + pj(n) , (A.4) 



p 3 {n + l) 

= 6tl(n)p 3 (n)q 3 (n) + 120t 3 (n)p 2 3 (n)[g 3 (n) + h 3 {n)\ 
+Ut 3 (n)q 3 (n)[g 3 (n) + h 3 (n)f + 88p 3 (n)[g 3 (n) + h 3 (n)} 3 
+4t 3 (n)p 3 (n)q 3 (n)[13g 3 (n) + 10h 3 (n)} + 78t 3 {n)p\{n) 
+Qp 2 3 (n)[A9g 2 3 (n) + 78g 3 {n)h 3 {n) + 29h 2 3 (n)} 
+2g 3 (n)[ll(73 3 (n) + 26 g\ (n) h 3 (n) + 19g 3 (n)h 2 3 (n) + 4h 3 3 (n)} 
+t 2 3 (n)q 2 3 (n) + 2t 3 (n)q 2 3 (n)\2g 3 (n) + h 3 (n)} + 26t 3 (n)p 2 3 (n)q 3 (n) 
+2p 3 (n)q 3 (n)[38g 2 3 (n) + 50g 3 (n)h 3 (n) + lhh 2 3 {n)\ 
+2p 3 {n) [llhg 3 {n) + 76h 3 (n)\ + 2t 3 (n)p 3 (n)q 2 3 (n) 
+q 2 3 (n)[Ag 2 3 (n) + Ag 3 (n)h 3 (n) + h 2 3 (n)] 
+2p 2 3 (n)q 3 (n)[31g 3 (n) + 18h 3 (n)} + A9p 4 3 (n) 

+2p 3 (n)q 2 3 (n)[2g 3 (n) + h 3 (n)\ + Up 3 3 (n)q 3 (n) + p 2 3 (n)q 2 3 (n) , (A.5) 
93 (n+ 1) 

= U4t 3 (n)p 3 (n)q 3 (n)[g 3 (n) + h 3 {n)] + 208t 3 {n)p 3 3 {n) 
+720p 2 3 (n)[g 3 (n) + h 3 {n)f + 56^(71) [<? 3 (n) + h 3 {n)f 
+24t 3 (n)q 2 3 (n)[g 3 (n) + /i 3 (n)] + 252t 3 (n)p 2 3 (n)q 3 (n) 
+2Ap 3 (n)q 3 (n)[25g 2 3 (n) + Mg 3 (n)h 3 (n) + I9h 2 3 {n)\ 
+104p 3 i (n)[17(? 3 (n) + 15/i 3 (n)] + 60t 3 (n)p 3 (n)q 2 3 (n) 
+24q 2 3 {n)[3g 2 3 {n) + 5<? 3 (n)/i 3 (n) + 2h\{n)\ 
+12p 2 3 (n)q 3 (n)[U0g 3 (n) + 89h 3 (n)\ + 972p 4 3 (n) + At 3 {n)q 3 {n) 
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+12p 3 {n)ql{n)[22g 3 {n) + I7h 3 {n)\ + 776pjj(n)g 3 (ra) 
+4ql(n)[4g 3 (n) + 3h 3 (n)} + 210p§(n)gf (ra) + 24p 3 (n)g|(n) + g|(n) . 



(A.6) 
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